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Preface 


The changes incorporated in this revised edition are contained 
solely in the table "Simulations of Sub-Functions (1729 SUBFNS)" 
on page 21. The table has been reorganized, some functions have 
been reformatted and simplified, known bugs have been corrected, 
and some new functions have been added. 


Robert A. Smith 
11 July 1975 


Preface to First Edition 


This set of very rough notes and tables was used as the working 
basis for the first Special Topics Seminar at Scientific Time 
Sharing Corporation. The Seminar was held 10 February 1975 in 
STSC's Washington, DC, Marketing Office. 


Some of the tables have been reorganized since their initial 
presentation at the Seminar. In particular, the table "Auxiliary 
Functions WA and PA" was completely rewritten, much to its 
betterment. 


The notes grew out of the author's deep and continuing interest 
in the subject and his desire to proselytize to others. 
Unfortunately, these notes are not self-contained. At the 
Seminar, much verbal accompaniment substituted for written prose 
to make the notes come alive. It is hoped that the reader will 
overcome this deficiency and make his own effort to bring life to 
them in the absence of an interactive seminar. 


While this work only scratches the surface in this field, already 
many problems can be solved with what is contained herein. The 
author believes that these ideas can be useful to virtually all 
APL programmers. Hopefully, some readers will be interested 
enough to further this work. 


Robert A. Smith 

Scientific Time Sharing Corporation 
Bethesda, MD 

17 February 1975 
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EQUIVALENCES FOR A@A®@B IN THE FORM OF A@B 
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(A<B)<A<C 
(A<B)<A<C 
(A<B)sA<C 
(A<B)2A<C 
(A<B)>A<C 
(A<B)#A<C 
(A<B)VA<C 
(A<B)AA<C 
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A¥B2C 
AVESC 
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(A#B)<A#C 
(A#B)SA2#C 
(A#B)=A2#C 
(A#B)2A2#C 
(A#B)>Axé 
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A<B<C 
AzB<C 
AAB<C 
A®B<C 
AsEsC 
A>BsC 
AVBsC 
A¥BSC 
ASB=C 
A>B=C 
AVB=C 
A¥B=C 
AsB2C 
A>B2c 
AVB2C 
A¥B2C 
A<B>C 
A2B>C 
AAB>C 
ARB>C 
A<B2#C 
A=BeC 
AAB#C 
AWwRZC 
A<BVC 
ASPVC 
A2PVC 
L£>PVC 
AVRVC 
AABVC 
A*BVC 
AWRVC 
A<BAC 
ASBAC 
A2BAC 
A>BAC 
AVBAC 
AABAC 
A¥BAC 
ANBAC 
A<B¥C 
ASB¥C 
A2BvC 
A>B¥C 
AVB¥*C 
AAB¥C 
ABYC 
A®BY¥C 
A<Bee 
ASBRC 
AZ2Bee 
A>Bel 
AVBeC 
AABRC 
A*BRC 
AWBHC 


(A<B)<A<C 


(A2B)sA2C 


(ASB )<ASC 
(A>B)<A>C 
(ASB)<SASC 
(A>B)<A>C 
(A<B)<A<C 
(A2B)<A2C 
(ASB)=A<C 
(A<B)#A<C 
(A<B)=A<C 
(A<B)#A<C 
(A>B)<sA>C 
(ASB)<ASC 
(A2B)sA2C 
(A<B)<A<C 
(A2B)<A2C 
(A<B)<SA<C 
(A>B)<A>C 
(ASB)sSASC 
(A<B)#A<C 
(A<B)=A<C 
(ASB )#ASC 
(ASB )=ASC 
(A<B)VA<C 
(ASB)VASC 
(A222 )AAzC 
(A>? )AAPC 
(AVE )VAVC 
(AAB)VAAC 
(A¥B)AA¥C 
(A®R)AARC 
(A<B)AA<C 
(ASB)AASC 
(A2P)vA2C 
(A>B)VA>C 
(AVB)AAVC 
(AAB)AAAC 
(A¥B)VA¥*C 
(AWB )VAAC 
(AB )AA¥C 
(AWB )AARC 
(AVB)VAVC 
(AAB)VAAC 
(A2R)AA2C 
(A>B)AA>C 
(A<B)VA<C 
(ASB)VASC 
(A¥B)VA¥C 
(AWB) VAKC 
(AVB)AAVC 
(AAB )AAAC 
(A2B)vA2C 
(A>B)VA>C 
(A<B)AA<C 
(ASB )AASC 
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(AVB)<AVC 
(A¥B)<A¥C 
(AAB)<AAC 
(AWB)<AKRC 
(AAB)SAAC 
(AXB)<AKC 
(AVB)sSAVC 
(A¥B)<A¥C 
(A>B)=A>C 
(A>B)#A>C 
(A2B)=A2C 
(A2B)#A2C 
(AWB) SANC 
(AAB)<AAC 
(A¥B)<A¥C 
(AVB)<AVC 
(A¥B)<A¥C 
(AVB)sAVC 
(AWB)<AKC 
(AAB)SAAC 
(A2B)#A2C 
(A2B)=A2C 
(A>B)#A>C 
(A>B)=A>C 
(A2B)wA2c 
(A>B)wA>C 
(A<B)¥A<C 
(A<B)¥A<C 
(A¥B)#A¥C 
(AWB )RARC 
(AVB)¥AVC 
(AAR)¥AAC 
(A2B)¥A2C 
(A>B)¥A>C 
(A<B)WA<C 
(ASR )WASC 
(AB )¥A¥C 
(AWB )¥AKC 
(AVB)WAVC 
(AAR )WAAC 
(AVB)¥AVC 
(AAB)¥AAC 
(AB )WA¥C 
(AWB )RARC 
(A<B)¥A<C 
(ASB)¥ASC 
(A>B)wA2C 
(A>B)wA>C 
(AVB)wWAVC 
(AAB)®AAC 
(A¥B )¥A¥C 
(A®B )¥ARC 
(A<B)WA<C 
(ASB)WASC 
(A2B)¥A2C 
(A>B)¥A>C 
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(A2B)>A2C 
(A<B)2A<C 
(A>B)>A>C 
(ASB)2ASC 
(A>B)2A>C 
(ASB)>ASC 
(A2B)2A2C 
(A<B)>A<C 
(AAB)=AAC 
(AAB)#AAC 
(AVB)=AVC 
(AVB)#AVC 
(ASB)2ASC 


(A>B)>A>C 


(A<B)2A<C 
(A2B)>A2C 
(A<B)>A<C 
(A2B)2A2C 
(ASB)>ASC 
(A>B)2A>C 
(AVB)#AVC 
(AVB)=AVC 
(AAB)#AAC 
(AAB)=AAC 
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(A¥B)>A¥C 
(AVB)2AVC 
(AWB)>AWC 
(AAB)2AAC 
(AWB )2AwWC 
(AAB)>AAC 
(A¥B)=A¥C 
(AVB)>AVC 
(AWB )=ANC 
(AWB )#AWNC 
(A¥B)=A¥C 
(A¥B)#A¥C 
(AAB)2AAC 
(AWB )>AnC 
(AVB)2AVC 
(A¥B)>A¥C 
(AVR)>AVC 
(A¥B)2A¥C 
(AAB)>AAC 
(AwB)2AWNC 
(A¥B)#A¥C 
(A¥B)=A¥C 
(AwB)#AWC 
(AWB)=AKNC 


Interpretations of all possible reductions of a 
Boolean vector where the result is 0 or 1 


Let B be any Boolean vector. 


REDUCTION 
i. A/B 
ae v/B 
3. w/B 
u, ¥/B 
5. </B 
6. </B 
‘- =/B 
8. 2/B 
9, >/B 

40. #/B 

ii. x/B 

S2« [/B 

ee L/B 

Th * /B 

15. |/B 

16. i/3 


INTERPRETATION 

0 if Oe€B; 1 otherwise. 

1 if 1e€B; 0 otherwise. 

>/B if RBe+N4i.N or WNpil;3 ~2|+/A\1=B | 
(the reverse parity of the number of leading 
1's in 8B) otherwise. 

>/B if Be+N=1N or WNp03 2/+/A\0=B 

(the parity of the number of leading 0's in 
B) otherwise. 

1 if B+e+N=.1N; 0 otherwise. 


0 if B+e+N4i1N3; 1 otherwise. 


~2|+/0=B (the reverse parity of the number 
of o's in 8B). 


~2|+/aA\0=B (the reverse parity of the number 
of leading 0's in B). 


2\|+/A\1=B (the parity of the number of 
leading 1's in 8B). 


2\|+/1=B (the parity of the number of 1's 
in 8B). 


See #1. 
See #2. 
See #1. 
See #8. 
See #5. 
See #6. 


'4+' means ‘is of the form’. 
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EIGHT DIFFERENT DEFINITIONS OF SCAN, WITH FXAMPLES 


R<« a\V #\15 -\15 
1. RCI] + a/ Iee+ve | 14 3 6 10 15 | *"4 2 
2. R(I] « a/ (7F-1)4+V | 15 1412 #9 5 | 3°72 4” 
3. P[I] « a/ (-F)4+v | 5 9 12 14 15 | 574 4” 
ub, RCI] + a/ (i-r)+vV | 1510 6 3 12 | ed a od 
5. P({I] + a/d Itv | 1 3 #6 10 15 | : ao on 
6. RCI) + a/Oo(7F-1)4+¥7 | 45 14412 #9 #5 | 3 2 4 
7, PCI] « a/d (-I)t+v | 5 9 12 14 15 | 5 41 4 
8. RCI] «+ a/o(i-ZT)+V | 1510 6 3 1 | 3 2 2 


IDFNTITIFS INVOLVING SCAN 


<\B «> ~s\~B 
<\B «> ~<\~B 
>\B +> ~>\~P <— =\v\B 
>\B «> ~2\~B <> #\A\B 


=\R <> ~z#\~B 
z\B +> ~=\~B 
vV\B «++ ~A\~B 
A\R +> ~vV\~B 
¥\B <> ~w\~B <> (2\R)=(V\P)S<\B 


#\B <> ~w\~B +> (>\P)#(A\B)<S\B 


MOF DF NF dD 
RPOOWde ® UI Ww 


Similarities Between Reduction and Scan 


of Boolean Vectors 


The table below illustrates notational similarities 
between APL expressions for reduction and scan of boolean 
vectors. The functions considered are all the primitive 
scalar dyadic relational and logical functions. Origin 
dependent expressions should be evaluated in origin l. 


Simpler statements exist for some of the expressions. 


For instance >/B «+ ~2|/+/A\~B ++ 2/811 and >/B ++ 
2|/+/A\B +> ~2|B10. The form using A\ was not chosen in 
order to point out the similarities with the expressions 
illustrating other reductions and scan. Also, A\ is one 


of the objects being illustrated and should not be used to 
define other expressions. The other equivalent form was not 
chosen since it has no simple extension to an expression for 
a scan uSing that function. 


REDUCTION SCAN 

A/B ++ (Br10)>pB A\B +> (B10)>19B 

v/B ++ (Bi1)s90B v\B <> (Bi11)<109B 

n/B ++ ((B10)<pB)=~2| w\B <«+ ((B10)<1pB)=~2| 
+/(B10)>19B +\(B10)>109B 

~/B <> ((Bi1)<pB)=2| w\B ++ ((B11)<1pB)=2\| 
+/(Bi1i)>19B +\(Bi1)>19B 

</B ++ (B11)=pB <\B +> (B11)=19B 

</B ++ (B810)#pB <\B ++ (B10)#19B 

=/B +> ~2|+/~B =\B ++ ~2|+\~B 

>/B +> ~2|+/(B11)>19B >\B ++ ~2/+\(B11)>19B 

>/B +> 2|+/(B10)>19B >\B «+ 2|+\(8B10)>10B 

z/B +> 2|+/B z\B +> 2/|+\B 


'++' means 'is equivalent to' 


Robert A. Smith 

Scientific Time Sharina Corp. 
7316 Wisconsin Avenue 
Bethesda, Md. 20014 
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ROTES 


R+v\B 
P+a\B 
R<+2\B 
R+>\B 
R+<\B 
R+s\B 
R+\P 
R+«w\F 


PAST ALGORITHMS FOR IMPLEMENTING RELATIONAL AND LOGICAL 


SCANS OTHER THAN #\ AND =\ 
: ALL COMPUTATION IN ORIGIN 13; 4 ++ (pB)-A 
A+ 14+Bii1 = R#(Ap0 0),d4pi 
A+ 1+Bi0 « R#(Api 1),4p0 
A+ 1+Bil = Re(Apd 1),4p~2]A 
A+ 1+P10 ™ Re(Api 0),4p 214A 
A+ 1+Pii = R+(Ap0 0),490 x IF AtoB THEN RCLA+1I)+1 
A+ 1+Bi0 = R+(Api 1),4p1 x IF A#toB TREN RCA+1J<+0 
A+ 1+Bii « R+(Ap0 1),4p 2/A = JF AtoB TREY RLA+1)+~2/A4 
A+ 1+Bi10 ® R+«(Api 0),4p~2/A = JF AtpB THEN RLA+1I)+ 2A 


ROBERT A. SMITH 
STSC 
SEPTEMBER, 1973 
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Function XVF: Expansion Vector for Boolean Right Arguments 


C1] 

C2] 

C3] 

C4] 

[5] 

C6] 

C7] 

C8] 

[9] 

[10] 
C11] 
C12] 
(i3] 
C14] 
C15] 
C16] 
C17] 
[18] 
C19] 
[20] 
C21] 
orn 
[23] 
[24] 
[25] 
[26] 
f27) 
[28] 
[29] 
[30] 
C31] 
(32] 
[33] 
C34] 
[35] 


VXVFCOIV 

V RL XVF B;C 
a EXPANSION YECTOR FOR BOOLEAN RIGHT ARGUMENT. 
a COMMENTS THRU [13] 
a CREATE AN EXPANSION VECTOR FOR A COORDINATE WITH LENGTH 
a EQUAL TO p<RIGHT ARGUMENT>. USE THE ELEMENTS OF THE LEFT 
a CYCLICALLY TO DETERMINE THE NUMBER OF ZEROES TO INCLUDE 
a NEXT TO ELEMENTS CORRESPONDING TO ONES IN THE RIGHT 
a ARGUMENT. NEGATIVE INSERTIONS GO TO THE LEFT, POSITIVE 
a TO THE RIGHT. 

& 

ee 

2g 

mR 

e 


ALGORITHM FROM DISCUSSION WITH BOB SMITH/STSC 
BRIAN HAGENBUCH STSC 30 OCT 74 


CLEAN UP THE ARGUMENTS. NOTE: NO ERROR CHECKING IS DONE. 
R+pB+,B * BB/ik 
L+(pP)pL 
a REMOVE LENGTH 0O INSERTIONS 
C+O4L & B+C/B 
[<+C/L 
C+0<L ma FLAGS ‘RIGHT' INSERTIONS 
a INJTIALIZE <R> TO 1 GREATER THAN THE LENGTH OF THE RESULT 
a T0 AVOID INDEX ERRORS. 
L+|L « Re(1+R++/L)pl 
a PUT O'S WHERE EACH INSERTION SHOULD BEGIN. 
B+B+C++\ 14+0,2L 
R{LBJ]+0 
a PUT O'S IN THE POSITIONS AFTER EACH INSERTION SHOULD END 
a EXCEPT WHERE THEY COINCIDE WITH BEGINNINGS, IN WHICH 
a CASE THE 1 SHOULD BE RESTORED. (SEE SCAN COMMENT BELOW. ) 
BeB+L 


RLB1]<+~R(B] 

a EQUALS SCAN TO CHANGE 1'S BETWEEN PAIRED O's TO-% . ces 

A (E.G. -..011011410010%11 1... BECOMES 
A 0001111011000 0...) 

a... AND DROP THE EXTRA POSITION FROM <R>. 

R+ 14+=\R 


V 


=F 


Sub-Functions in 4PL 


The Problem 


The coordinates of matrices and higher dimensional arrays 
serve to delimit the values of the array into distinct “lines” 
(e.g., the rows of a matrix). Monadic APL primitives (e.g., 
reduction and scan) then may be applied independently to each 
distinct line to produce various results. 


Occasionally, however, it becomes necessary to apply a monadic 
function independently to successive parts of a vector (called the 
argument vector). These parts are analogous to the lines of a 
higher dimensional array. For example, the problem may be to sum 
reduce distinct parts of a numeric vector. If the argument vector 
were 110 and if it were broken into several parts as follows: 


4258 My > 6,4,8 9 10 


then the desired result would be 6 4 26 19 . Previously, this 
kind of problem has been solved by first expanding and reshaping 
the argument vector into a matrix. This transformation allows the 
programmer to take advantage of the delimiting property of the 
coordinates of the matrix and apply the monadic function to the 
entire matrix to obtain the result. If the argument vector were 
stored in B and the lengths of the successive parts stored in 
A, a solution is +/((pA).[/A)p(,A°.21f /A)\B . 


As a general method, the argument vector-to-matrix approach 
has several drawbacks. First, the expansion result may not fit in 
the workspace. Since the length of the expansion result is the 
number of parts times the length of the longest part, it only takes 
one long part and many small parts to create a nuisance. WS FULL 
is my verv least favorite error message. Second, the expansion 
process may "pad" each part (on the way to creating a row of the 
matrix) with the "wrong" value. The zeros which expansion uses as 
fillers for numeric structures won't bother +/ but *x</ and A/ 
won't like it. The programmer can get around it, but not without 
expending rather more CPU seconds (the code is easy, though). 
Third, many CPU seconds can be used by this approach due to all the 
fill elements inserted by the expansion. Fourth and last, the 
monadic function may not be defined for matrices. While I have not 
seen the need to perform a grade up on successive parts of a 
vector, it could happen. 


The following sections discuss an alternate approach to this 
kind of problem. The method consists of characterizing the concept 
of "parts" and then combining that with the monadic function and 
the argument vector to produce the result. 


= 


The “sub"-operator 


A simple and unique characterization of any partitioning of 
an argument vector is as a Boolean vector of the same length, with 
each 1 corresponding to the beginning of each part and the 


following 0's corresponding to the remaining elements in that 
part. 


The partition of 110 used above would be characterized as: 


3 4 > 6 9 10 
0 1 1 0 


i. 2 7 8 
1 0 0 0 i 0 

It is easy to see that this characterization is unique and 
that every Boolean vector of the same length as the argument 
vector and whose first element is 1, represents a partition of 


the argument vector. There are 2*N-1 such partitions of length 
N. 


The "sub"-operator takes a monadic function and produces a 
dyadic partitioning function. The resulting function is dyadic 
because it now needs a Boolean vector (left argument) to specify 
the partition. 


For simplicity, a new symbol is used for this operator. It 
is not the goal of this paper to suggest that this operator or 


symbol be implemented; it is introduced only to simplify the 
presentation. 


Definition: sub-operator + (which Jim Ryan calls the 
—— "dagger" symbol) 


Let a be any monadic function, and A and B any 
equal length vectors where A is a Boolean vector whose 


first element is i. 
R+A+aB 


is the result of applying a to each part of B where A 
defines the partition. Specifically, 


R+10 
FOR JIeit/A DO R+R,a(#\A\(1+/A )eI+ 0 4)¢B%4 
The shape of R is the sum of the number of elements in the 
result of applying «a to each part of B. If a is defined for 


higher dimensional arrays, then so is ¢a and in the same way. 
Notationally, this would be written as R+AtalXK] B 


The above example of a partition and sum reduction of 110 
would be written as 


1001410001 0 t+/110 


or with B+110 and A+3 14 2 = (the partition lengths) as 


(CipB)eitt+\ 140,A)t+/B. 
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An Application 


The need which motivated this entire effort was to write a 
function which produced a cross reference of all identifiers used 
in an arbitrary function. To do this, all elements within quotes 
or comments must be identified. First, a listing of the function 
was obtained which duplicated the display form (VFNNAME[U]V) of 
the function as a character vector, with CR's separating the 
lines of the function, CR's in character constants followed by 
six blanks, etc. A primitive function (AFD) exists on the 
APL*PLUS System to produce this character vector, given the name 
of the function. Next, the elements of the function which were 
contained between quotes had to be identified. Here is where the 
problem began. A change to the APL*PLUS System had eliminated the 
annoying practice of the system doubling all quotes in a comment 
line. This change had made it entirely possible for a function 
line to contain an odd number of quotes, @€.9., 


C3] a WHAT'S UP DOC? 


Prior to this change, a simple #\V='''' would have sufficed. 
Not only that, but another change to our system had begun the 
useful practice of allowing comments to appear as the rightmost 
statement, even on lines with executable statements, e€.9., 


C3] [+0 « LIM*«(FLIM 10)(2]-1 a SET INDEX VAR AND LIMIT. 


One could no longer rely upon a comment to appear at a fixed 
position in a function line. So, a comment symbol could appear in 
quotes; an unmatched quote could appear in a comment; and a 
comment need not be the only statement on a line. The problem of 
identifying all elements within quotes or comments looked sticky. 


However, a working solution was found which later was 
recognized to be equivalent to: 


FL+1 AFD ‘FNNAME' a VFNNAME(OJV FUNCTION LISTING 
DEL«(FL='"(')A 140,FL=CR ® DEL(11J+1 a PARTITIONING VECTOR 
Q+DELt#\FL='''* & RES*QVDELtV\Q<FL='a' 


Since the vector display form of the function (the result of 
1 AFD ) follows each CR in a character constant with six blanks, 
the i's in DEL correspond to the left bracket which begins 
each line of the function. This partitioning vector in 
combination with the sub-operator allows one to work independently 
on each line of the function listing. In particular, Q@ isa 
Boolean vector which selects all elements (even if in a comment) 
which are to the right of an odd number of quotes in a function 
line. The 1's in Q<FL='a' select the lamp symbols that are 
not within quotes, and DEL+v\ of this selects all elements 
within comments. The final result selects all elements within 
quotes or comments. 


 & 


Tools and Techniques 


All this discussion would be only academic if there were no 
efficient (read “non-looping") implementation of these ideas. 


In all, the set of monadic primitives which do something 
useful with a vector argument are the 13 primitive scalar 
monadic functions, reduction, scan, grade up, grade down, and 
reversal. For the 13 PSMF's, AtaB ++ aB}; if the function 
argument involves no interaction between the elements, the 
partition has no effect. Of the 21 PSDF's, 7 of the more 
commonly used functions have been simulated and are listed in the 
Appendix in both scan and reduction form. Sub-functions which 
simulate grade up, grade down, and reversal are also included more 
for completeness than utility. In most cases, each function 
requires only about the same amount of storage in which to execute 
as does the function when not used in the context as a 
sub-function. Each function is non-looping and works only on 
vector arguments. All functions work in either origin with +4 
and +¥ returning results dependent upon the origin. 


A tool used in implementing these functions is the function 
NA (cf. [1]). The function works as follows: 


Ta! VA Vo «++ Va 1+(a/10),V . 


The left argument must be a character scalar symbol for a 
primitive scalar dyadic function with an identity element. 
Essentially, WA compares V with itself shifted one to the 
right (negatively, hence the W in WA), using the function 
symbol a. The shift of V drops off the last element and 
catenates the identity element of a to the beginning. For 
example, ‘'#' WA V ++ Vz 1+0,Y. A companion function, PA , 
shifts to the left (positively) and fills the hole (this time at 
the end) in the shifted vector with again the identity element of 
ae ‘at PA V ++ VaitV,a/10 . 


In the abstract, WA and PA are operators which take a 
primitive scalar dyadic function with identity element and return 
a monadic function which performs a "differencing" (hence the A 
in their names). WA is similar to the function P4JR in 
reference [1], where Halpern also notes that that function is the 
inverse of certain scans. For example, ‘'+' WA and #\ are left 
and right inverses of each other as are ‘=' WA and =\ ; ‘'-' WA 
and +\ 3; and ‘'#' WA and x\ . 


The techniques used in implementing the sub-functions are 
interesting enough to warrant stepping through some examples. 
Note that for a partitioning vector A , the i's in A_ are the 
starting points of each subdivision, while the i's in 1A are 
the corresponding endpoints. 
is A¢t+\B ++ +\B-A\'-' NWA A/t+\ 140,B8 . 


Let 


~1h= 


A 1 0 9) 1 1 0 0 0 1 0 
B 1 2 3 m 5 6 7 8 9 10 
14. +\ 140,B 0 1 3 6 10 15 21 28 36 45 
2. af 0 6 10 36 
3. ‘'-' WA 0 6 4 26 
uy, 6 6OA\ 0 0 0 6 4 a) 0 0 26 0 
5. B- 1 2 3 ~ 2 1 6 7 8 “SPF EO 
6. +\ | 3 6 4 - Sees op items 3. Ube s 9 19 


Steps 1 and 2 compute, for the four subdivisions of 8B , the 
cumulative sum of all the previous subdivisions. Step 3 
isolates, for each subdivision, the contributions of the 
immediately preceding subdivision. Step 4 puts the result of 
Step 3 into the positions of the beginning of each subdivision, 
and Step 5 subtracts this from 8B. The result at this point is 
similar to B but with the initial value in each subdivision 
short by the sum reduction of the immediately preceding 
subdivision. When Step 6 is executed, each subdivision will 
have the bias from any previous subdivision accounted for 
correctly. 


Pleasantly, At+#\B falls out of A++\B because for a Boolean 
vector B , *\B ++ 2/+\B . Thus, 


Atze\B <> 2/At+\B 

++ 2/+\B-A\'-"' WA A/+\ 140,B 
z\2|B-A\'-"' NA A/+\ 14+0,B 
4\Be2|A\'-' WA A/+\ 14+0,B 
4\BeA\2|'-"' HA A/+\ 140,B 
x\BeA\'s' HA 2|A/+\ 14+0,B 
z\BeA\'tz" WA A/2|+\ 14+0,B 
z\BeA\'z' WA A/#\ 140,B . 


pi ddys 


A+=\B is derived from the identity =\B +> ~#\~B . 
II. At+/B +> '-' WA(1OA)/+\B © 


This code finds the cumulative sum at the endpoints of each 
subdivision and pairwise differences these values to obtain the 
result. The Appendix contains another function to do A*¢+/B 
which requires much less workspace than the above method when B 
also is a Boolean vector. 


III. Atv\B ++ #\(AVB)\'#' NAC(AVB)/B . 


Let 

A 100 1 1000 1 0 

B 010 0 1001 0 0 
1. <AVB 11 0 1 1001 1 0 
as (AVB)/B Oo 1 0 1 1 0 
- 's' WA 01 1 1 0 1 
4, (AVB)\ 01 0 1 1000 1 0 
5. #\ O11 0 5 ee ae ae | 0 0 


Step 1 points out that the important information in A and B 


wt es 


is where either of them is 1. Where both of them are 0 , the 
value in the corresponding position in the result will be equal to 
the value in the immediately preceding position in the result. 

The leading i's in each set of consecutive i's in Step 2 
indicate a 1 in 8B where the preceding elements in that same 
subdivision of B were all 0. This position in B represents 
a point at which we wish to begin smearing i's to at least the | 
end of that subdivision. The actual point at which we should stop 
smearing i's is where the initial element in some subdivision of 
B is 0. These stopping points are exactly the points in 8B 
corresponding to the 0's in Step 2. Step 3 results ina 1 
at both the start and stop points of each smear of i's. Step 4 
puts these values back into place relative to 8B and Step 5 
performs the smear. 


The relationships between the various steps in this solution 
are quite intriguing. Steps 2 and 4 are a compression and 
expansion using the same left argument. Steps 3 and 5 are 
inverses of each other. 


A+a\B is derived using the identity A\B ++ ~V\~B . 


Reference: 
[1] Halpern, Michael M. "Studies in APL: Algebra, Scan, 
Arithmetic, Permutations" IBM Philadelphia Scientific 
Center, Tech Report No. 320-3023, June 1973. 
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ORIGIN CLASSIFICATIONS FOR PRIMITIVE FUNCTIONS 
ORIGIN SENSITIVE: []1?7498& 


ORIGIN FREE: ALL OTHERS. 


ORIGIN CLASSIFICATIONS FOR DEFINED FUNCTIONS 


ORIGIN FREE: THE PROBLEM WARICH THE FUNCTION SOLVES HAS TRE SAME 
ANSWER INDEPENDENT OF THE SETTING OF THE ORIGIN, AS DOES TRE 
FUNCTION. TRE FUNCTION USES ONLY ORIGIN FREE PRIMITIVES AND 
WORKS CORRECTLY IN EITHER ORIGIN. 


ORIGIN SENSITIVE: THE PROBLEM WHICH THE FUNCTION SOLVES HAS A 
DIFFERENT ANSWER DEPENDING UPON THE SETTING OF TRE ORIGIN, AS 
DOES THE FUNCTION. THF FUNCTION USES ORIGIN SENSITIVE PRIMITIVES 
BUT WORKS CORRECTLY IN EITHER ORIGIN. | 


ORIGIN INDEPENDENT: THE PROBLEM WHICH THE FUNCTION SOLVES HAS THE 
SAME ANSWEP INDEPENDENT OF THE SETTING OF TRE ORIGIN, AS DOES THE 
FUNCTION. TRE FUNCTION USES ORIGIN SENSITIVE PRIMITIVES BUT 
WORKS CORRECTLY IN EITHER ORIGIN. 


ORIGIN 0 (OR 1) DEPENDENT: THE PROBLEM WHICH THE FUNCTION SOLVES MAY 
OR MAY NOT DEPEND UPON THE SETTING OF THE ORIGIN. THE FUNCTION 
USES ORIGIN SENSITIVE PRIMITIVES AND WORKS CORRECTLY ONLY IN THE 
SPECIFIED ORIGIN. 


ROBERT A. SMITH 
STSC 
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